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We study theoretically the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction in one- and two- 
dimensions in presence of a Rashba spin-orbit (SO) coupling. We show that rotation of the spin of 
conduction electrons due to SO coupling causes a twisted RKKY interaction between localized spins 
which consists of three different terms: Heisenberg, Dzyaloshinsky-Moriya, and Ising interactions. 
We also show that the effective spin Hamiltonian reduces to the usual RKKY interaction Hamiltonian 
in the twisted spin space where the spin quantization axis of one localized spin is rotated. 



There has been a great deal of interest in the field of 
spintronics where spin degrees of freedom of electrons 
are manipulated to produce a desirable outcome^'^. Em- 
inent examples are given by the giant magnetoresistance 
(GMR) effect'^''*'^ and the interlayer exchange coupling in 
magnetic multilayer o^i^i^ . The interlayer exchange cou- 
pling is explained in the context of Ruderman-Kittel- 
Kasuya-Yosida (RKKY) interaction^' equivalently, or 
in terms of spin-dependent electron confinemen1iii*i^. 
The RKKY interaction is an indirect exchange interac- 
tion between two localized spins via the spin polarization 
of conduction elect ron3iii^*iS*i&. 

Recently, much attention has been focused on the 
effect of the Rashba spin-orbit (SO) coupling in two- 
dimensional electron gases (2DEG)^^. Investigation of 
the Rashba effect of 2DEG in semiconductor heterostruc- 
tures has been stimulated by the proposition of a spin 
field effect transistoui^. It has been established the 
Rashba SO coupling can be controlled by means of a 
gate voltage^^-^"-^^-^^. The Rashba effect has also been 
observed in 2DEG formed from surface states electrons at 
metal surfaces such as Au(lll)SiSia5i2SiSl, Li/W(110) or 
Li/Mo(110)^. It has also been found that confinement 
of the surface state due to atomic steps on vicinal sur- 
faces leads to quasi one-dimensional (ID) surface states, 
which also exhibit the Rashba effecliSSiSSi^. 

Usually the RKKY interaction yields a parallel or an- 
tiparallel coupling of localized spins (Heisenberg cou- 
pling). However, if spin of conduction electrons pre- 
cesses due to the spin-orbit coupling, it can be possible 
to produce a non-coUinear Dzyaloshinsky-Moriya (DM) 
coupling of localized spins^SiSi^. In this paper, we in- 
vestigate the RKKY coupling between localized spins em- 
bedded in a ID- or 2DEG with Rashba SO couphng. We 
show that rotation of the spin of conduction electrons due 
to the Rashba SO coupling causes a twisted RKKY in- 
teraction between localized spins which consists of three 
different terms: Heisenberg, Dzyaloshinsky-Moriya, and 
Ising interactions. We point out that a perturbative 
treatment of the SO coupling as is usually done^^-^"^-^^ is 
valid only for small distances between the localized spins; 
in this case the DM and Ising terms are respectively linear 



and quadratic with respect to the SO coupling strength. 
In the limit of large distances, a non-perturbative treat- 
ment of the SO coupling is necessary, and one obtains 
DM and Ising terms that have the same oscillation am- 
plitude as the Heisenberg term, independently of the SO 
coupling strength. This peculiar behavior of the twisted 
RKKY interaction for a pair of localized spins can be 
explained by introducing a twisted spin space where the 
spin quantization axis of one of the localized spins is ro- 
tated. 

We consider a system consisting of two localized spins 
embedded in a ID- or 2DEG with a Rashba type spin- 
orbit couplingii. The Hamiltonian for the conduction 
electrons is given by 

Hn = - — V^ + a{-ihV y. z) ■ (T, (1) 
2m 

where a represents the strength of the spin-orbit cou- 
pling, i is a unit vector along the 2;-axis, and cr is the 
vector of Pauli spin matrices. We assume that the con- 
duction electrons are confined in a wire along the x-axis 
(one-dimensional system) or in the x — y plane (two- 
dimensional system). The direction of the effective elec- 
tric field of spin-orbit coupling is taken to be along the 
z-axis for both one- and two-dimensional systems. 

Since the Hamiltonian Hq commutes with the momen- 
tum operator ~ih^ , the wave vector A; is a good quantum 
number. The Green function of the conduction electrons 
in the real space can be expressed as 



G {R; z) = j d^k e^'^-^ G(fc; z) 



(2) 



where D = 1 or 2 is the dimension of the system and the 
Green function in the momentum space is given by 



G{k-z) 



21,2 



2m 



■ctq + a {k X z) ■ cr 



(3) 



Here (To is the (2 x 2) unit matrix in the spin space of 
conduction electrons. 

The localized spins are denoted by Si and S2 and lo- 
cated at positions Ri and R2, respectively. The cou- 
pling between conduction electrons and localized spins is 
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expressed as the s-d interaction Hamiltonian 
Hi = J Sir- R,)S, ■ a, 



defined as 



(4) 



i=l,2 



where J represents the strength of the s-d interac- 
tion. Note that J has the following dimensionality: 
(energy) x (length) . 

The total Hamiltonian is given by the sum of Hq and 
Hi. We assume that the coupling constant J is so small 
that we can treat Hi as a perturbation on Hq. The 
RKKY interaction between Si and S2 is calculated from 
the second order perturbation theory as 



rrRKKY 
-"1,2 



--Im /" deTT[{Si-(T)G{Ri2;e + iO+) 
(S2•c^)G(-i^l2;£ + ^0+)] , 

(5) 

where is the Fermi energy, R12 = Ri — R2, iO^ rep- 
resents an infinitesimal imaginary energy, and Tr means 
a trace over the spin degrees of freedom of conduction 
electronsiSi^^. 

Let us first consider the one-dimensional case. The 
Green function of conduction electrons in the momentum 
space is expressed as 



G{k;e + iO+) 



21,2 



~2m 



■do ^ akcr. 



(6) 



After some algebras, Eq. © can be written as 

G{k;e + iO+) ^ Go{k;e)<Jo + Gi{k;e)ay, (7) 
where the diagonal and off-diagonal Green functions are 



Goik;e) 



1 



h"^ k^-k^ - 2kkR + ^0^ 

1 

^k^-k^ + 2kkR + iO'^ 



Tfl 1 

Gi{k-e) = [kl-k^ -2kkR + i{)+ 

1 

~ /c| - P + 2fcfci? + iO+ 

with fc^ = 2me/h^ and kn = ma /It?. 

A straight-forward contour calculation gives 

G{R-e + i{)+) = Go{R;e)cro + Gi{R;e)ay, 

where 

GoiR;e) = -i- 

Gi{R;e) 



e"?l^lcos(/fcfl,i?). 



™ -e'«l^lsin(fcfli?) 



ft2(g + iO+) 



with 



2m 2 



(8) 



(9) 



(10) 

(11) 
(12) 

(13) 



Substituting Eqs. (CHll-OSl into Eq. 0, and using 
the relations Go(— i?;e) — Go{R;e) and Gi(— i?;e) — 
~Gi{R;e), we have 



-"1,2 



RKKY ^ _±j2j^^^ 



Tr{(Si-cr)(52-'T)} / ^ deGoiRi2;sf + Ti {{Si -a) ay [82 -(t)} f \e Gi{Ri2;e)Go{Ri2;e) 

J —00 J —00 

Tr{(5i-cr)(52-<T)aJ / ^ de Go(i?i2; £)Gi(i?i2; e) - Tr {(Si • *t) CTj, (^2 • <t) aj / ^ de Gi{Ri2\ef 

J —00 J —00 



I 

The traces over the spin operators can be calculated by Thus, using Eqs. Hll|l -H18 |l we find 
using the relation {A-cr){B-cr) = {A-B)ao+i{AxB)-cr, 
repeatedl}i2&: 



H^^^^" = Fii\Ri2\) 



COs{2kB,Rl2)Sl-S2 



sin(2fcfli?i2)(S'i X S2)y + {1 - cos{2kRRi2)} SfS^ 



Ty{{Si-ct){S2-(t)} = 2Si-S2, (15) 

Tr {(Si • cr) ay (82 ■ cr)} = -2i{Si x S2)y, (16) 

Tt{{Si- <t) {S2 ■(T)ay}^ 2i(Si X S2)y, (17) 

Tr{{Si ■ a) ay {S2 ■ (t) ay} = 2(25f^| - Si ■ S2). (18) 



where the range function i^i(|i?|) is defined as 

t2 2 f^F g^2iq\R\ 



{q + ?;o^ 



(19) 

2 de. (20) 
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Performing the change of variable e ^ q and using stan- 
dard complex-plane integration techniques, one eventu- 
ally obtains: 



Si{2qF\R\) 



(21) 



where qp = y^77iepyft^ + fc|j and Si( ) is the sine in- 
tegral function"^^. The range function of Eq. H21(l is 
the same form as that of the usual one-dimensional 
RKKY interactioE—i^ except that the Fermi wave vector 



fc_F(= ^/2rnepjh?) is replaced hy qp . 

As shown in Eq. H19|) . the resulting RKKY interac- 
tion consists of three physically quite different interac- 
tions: Heisenberg, Dzyaloshinsky-Moriya and, Ising in- 
teractions. The Heisenberg and Ising couplings favor 
a collinear alignment of localized spins. On the con- 
trary, the DM coupling favors a non-coUinear alignment 
of localized spins. For distances (more precisely, for 
kB.\Ri2\ ^ 1), the DM and Ising terms are, respectively, 
linear and quadratic in the Rashba SO coupling a; this 
corresponds to the result obtained from a perturbative 
treatment of the SO coupling. However, for large dis- 
tances (fcii;|i?i2| ^ 1), a perturbative treatment of the SO 
coupling would completely fail: indeed, in this regime, 
one finds that the DM and Ising couplings oscillate with 
the same amplitude as the Heisenberg term. 

This peculiar twisted coupling of localized spins can be 
easily understood by introducing the twisted spin space 
where the spin quantization axis of the second localized 
spin S2 is rotated by an angle 612 = 2fcfli?i2 around the 
j/-axis. The spin operators for the second localized spin 
in the twisted spin space are given by 



SUO12) = cos 5^ + sin 012 S^, 



(22) 
(23) 

(012 ) = cos 012 5| - sin 012 . (24) 



51(012) = SI 



From the above equations, one can easily show that the 
inner product of Si and ^2(012) is 

Si ■ ^2(012) = COS 012 Si -52 + Sin 012 (5i X S2)y 

-f {l-COS0l2}S'fS'2^ ^ ' 

SO that the RKKY interaction of Eq. (|19|) can be ex- 
pressed as 



H 



RKKY 
1,2 



Fi(|i?i2|)Si -52(012). 



(26) 



Eq. (|26|l shows that in presence of spin-orbit coupling the 
RKKY interaction results in a collinear coupling of local- 
ized spins in the 0i2-twisted spin space. 

Next we address the two-dimensional case. The Green 
function of conduction electrons now takes the form 

G{R;e + iO+) = Go{R;e)ao + Gi{R;e){z xR)-(t, (27) 

where R = \\R\\ and R = R/R is the unit vector parallel 
to R. The Green functions Go(i?;e) and Gi{R;e) can 



be calculated in the similar way as in Ref.4Q. A straight 
forward calculation yields 



GoiR;e)-- — 



Gi{R;s) 



im 




+ kR 


+ iO+)R 


+ ( 




kR + 


iO+)R] 


m 


(i+t) 


+ kR 


+ iO+)R 


-( 


:-^s)<'[(,- 


kR + 


iO+)R] 



(28) 



(29) 



where i?o^''[ ] and H^^\ ] are Hankel functions^. Here- 
after, we restrict ourselves to the region qR ^ 1 and 
kR <C q- In this case we can use the asymptotic form of 
Hankel's functions^, 

J?W(z):^,/Ae<--^-f) (|z|->oo). (30) 

V TTZ 



(«^-^)cos(fcfli?), (31) 



Thus we have 

GdR;e) = ^^±=e'("'-i) sin (k„R) . (32) 

Without restriction, we can take the coordinate system 
so that the vector R is aligned with the a;-axis, i. e., 
R = R X. Then the Green function takes the form 



G{R; e + iO+) = Go{R; e)ao + Gi{R; e)ay 



(33) 



The RKKY interaction can be obtained in the similar 
way as one-dimensional system, and one gets (for qpR^ 
1) 



^RKKY ^ F2{Ri2)Si ■ ^2(012), 

where the range function F2{R) is given by 

,P m sin {2qFR) 



F2{R) 



27r2 ft2 



(34) 



(35) 



Eq. (|34|l is the same as the usual two-dimensional RKKY 
interaction'^-'^- except that kp and 5^2 are replaced by 
qp and 52(012), respectively. It is reasonable that the 
twisted coupling of two localized spins takes the same 
form 51-52(012) as for the ID system, because for qpR ^ 
1 a scattering wave of 2D system behaves like a plane 
wave. 

The 2DEG from surface states at metallic surfaces are 
good candidates for investigating the twisted RKKY in- 
teraction. Experiments on Au(lll) surfaces have yielded 
qp ~ 0.17 and kR ~ 0.012 A"! ^; in complete 
analogy to the effect of negative gate voltage for 2DEGs 
at a semiconductor heterojunctioniSiSS, the adsorption 
of a noble gas (e.g., Xe) produces an effective repul- 
sive potential (because of the Pauli exclusion principle) 



4 



and leads to a decrease of the average Fermi wavevector 
{qf ^ 0.155 A-i) and an increase of the Rashba spUt- 
ting {kji ~ 0.015 A~^)2&. Furthermore, quasi- ID sur- 
face states can also be obtained for vicinal surfaces, such 
as Au(788) and Au(23 23 21), with comparable Rashba 
splittingS2i3fli2ii. The distance between magnetic adatoms 
deposited on such surfaces can be controlled either di- 
rectly by atom manipulation using the tip of a scanning 
tunneling microscope, or by exploiting self-organization 
processes. Such systems therefore constitute a versatile 
laboratory to investigate surface states mediated RKKY 
interactions under the influence of the Rashba effect. 
With kn ^ 0.015 A'^ and a distance i? = 10 A the 
twist angle 9 — 2kRR is of the order of 17°, which is 
quite sizeable. In particular, due to the twisted nature of 
the RKKY interaction, very interesting frustration phe- 
nomena may be anticipated. 

In semiconductor heterostructures the twisted RKKY 
interaction may also be of great interest, in partic- 
ular in view of possibilities for manipulating entan- 
glement between spins of quantum dots connected by 
a wire with Rashba SO couplingi^, as needed for a 



spin-based solid-state quantum computeriaaiaa. Consid- 
ering the value of the Rashba coupling reported for 
Ino.53Gao.47As/InQ.52Alo.48As heterostructureiS one ob- 
tains that the twist angle 9 can be controlled from 9 = tt 
to 61 = 37r/2 by a gate voltage for R — 400nm. 

In conclusion, we have studied the twisted RKKY 
interaction in one- and two-dimensions in presence of 
Rashba spin-orbit coupling. We have also shown that 
in the twisted spin space where the spin quantization 
axis of one localized spin is rotated, the twisted RKKY 
interaction is expressed in the same form as the usual 
RKKY interaction. The angle 9 between the localized 
spins can be controlled by the distance between localized 
spins R and/or the strength of the spin-orbit coupling of 
conduction electrons. 
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